Bayesian statistical analysis of bacterial diversity

Bacteria play an important role in many ecological systems. The molecular characterization of bacteria using either cultivation-dependent or cultivation-independent methods reveals the large scale of bacterial diversity in natural communities, and the vastness of subpopulations within a species or genus. Understanding how bacterial diversity varies across different environments and also within populations should provide insights into many important questions of bacterial evolution and population dynamics. This thesis presents novel statistical methods for analyzing bacterial diversity using widely employed molecular fingerprinting techniques. The first objective of this thesis was to develop Bayesian clustering models to identify bacterial population structures. Bacterial isolates were identified using multilous sequence typing (MLST), and Bayesian clustering models were used to explore the evolutionary relationships among isolates. Our method involves the inference of genetic population structures via an unsupervised clustering framework where the dependence between loci is represented using graphical models. The population dynamics that generate such a population stratification were investigated using a stochastic model, in which homologous recombination between subpopulations can be quantified within a gene flow network. The second part of the thesis focuses on cluster analysis of community compositional data produced by two different cultivation-independent analyses: terminal restriction fragment length polymorphism (T-RFLP) analysis, and fatty acid methyl ester (FAME) analysis. The cluster analysis aims to group bacterial communities that are similar in composition, which is an important step for understanding the overall influences of environmental and ecological perturbations on bacterial diversity. A common feature of T-RFLP and FAME data is zero-inflation, which indicates that the observation of a zero value is much more frequent than would be expected, for example, from a Poisson distribution in the discrete case, or a Gaussian distribution in the continuous case. We provided two strategies for modeling zero-inflation in the clustering framework, which were validated by both synthetic and empirical complex data sets. We show in the thesis that our model that takes into account dependencies between loci in MLST data can produce better clustering results than those methods which assume independent loci. Furthermore, computer algorithms that are efficient in analyzing large scale data were adopted for meeting the increasing computational need. Our method that detects homologous recombination in subpopulations may provide a theoretical criterion for defining bacterial species. The clustering of bacterial community data include T-RFLP and FAME provides an initial effort for discovering the evolutionary dynamics that structure and maintain bacterial diversity in the natural environment.

Methicillin-resistant Staphylococcus aureus in Finland: recent changes in the epidemiology, long-term facility aspects, and phenotypic and molecular detection of isolates

Staphylococcus aureus is one of the most important bacteria that cause disease in humans, and methicillin-resistant S. aureus (MRSA) has become the most commonly identified antibiotic-resistant pathogen in many parts of the world. MRSA rates have been stable for many years in the Nordic countries and the Netherlands with a low MRSA prevalence in Europe, but in the recent decades, MRSA rates have increased in those low-prevalence countries as well. MRSA has been established as a major hospital pathogen, but has also been found increasingly in long-term facilities (LTF) and in communities of persons with no connections to the health-care setting. In Finland, the annual number of MRSA isolates reported to the National Infectious Disease Register (NIDR) has constantly increased, especially outside the Helsinki metropolitan area. Molecular typing has revealed numerous outbreak strains of MRSA, some of which have previously been associated with community acquisition. In this work, data on MRSA cases notified to the NIDR and on MRSA strain types identified with pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and staphylococcal cassette chromosome mec (SCCmec) typing at the National Reference Laboratory (NRL) in Finland from 1997 to 2004 were analyzed. An increasing trend in MRSA incidence in Finland from 1997 to 2004 was shown. In addition, non-multi-drug resistant (NMDR) MRSA isolates, especially those resistant only to methicillin/oxacillin, showed an emerging trend. The predominant MRSA strains changed over time and place, but two internationally spread epidemic strains of MRSA, FIN-16 and FIN-21, were related to the increase detected most recently. Those strains were also one cause of the strikingly increasing invasive MRSA findings. The rise of MRSA strains with SCCmec types IV or V, possible community-acquired MRSA was also detected. With questionnaires, the diagnostic methods used for MRSA identification in Finnish microbiology laboratories and the number of MRSA screening specimens studied were reviewed. Surveys, which focused on the MRSA situation in long-term facilities in 2001 and on the background information of MRSA-positive persons in 2001-2003, were also carried out. The rates of MRSA and screening practices varied widely across geographic regions. Part of the NMDR MRSA strains could remain undetected in some laboratories because of insufficient diagnostic techniques used. The increasing proportion of elderly population carrying MRSA suggests that MRSA is an emerging problem in Finnish long-term facilities. Among the patients, 50% of the specimens were taken on a clinical basis, 43% on a screening basis after exposure to MRSA, 3% on a screening basis because of hospital contact abroad, and 4% for other reasons. In response to an outbreak of MRSA possessing a new genotype that occurred in a health care ward and in an associated nursing home of a small municipality in Northern Finland in autumn 2003, a point-prevalence survey was performed six months later. In the same study, the molecular epidemiology of MRSA and methicillin-sensitive S. aureus (MSSA) strains were also assessed, the results to the national strain collection compared, and the difficulties of MRSA screening with low-level oxacillin-resistant isolates encountered. The original MRSA outbreak in LTF, which consisted of isolates possessing a nationally new PFGE profile (FIN-22) and internationally rare MLST type (ST-27), was confined. Another previously unrecognized MRSA strain was found with additional screening, possibly indicating that current routine MRSA screening methods may be insufficiently sensitive for strains possessing low-level oxacillin resistance. Most of the MSSA strains found were genotypically related to the epidemic MRSA strains, but only a few of them had received the SCCmec element, and all those strains possessed the new SCCmec type V. In the second largest nursing home in Finland, the colonization of S. aureus and MRSA, and the role of screening sites along with broth enrichment culture on the sensitivity to detect S. aureus were studied. Combining the use of enrichment broth and perineal swabbing, in addition to nostrils and skin lesions swabbing, may be an alternative for throat swabs in the nursing home setting, especially when residents are uncooperative. Finally, in order to evaluate adequate phenotypic and genotypic methods needed for reliable laboratory diagnostics of MRSA, oxacillin disk diffusion and MIC tests to the cefoxitin disk diffusion method at both +35°C and +30°C, both with or without an addition of sodium chloride (NaCl) to the Müller Hinton test medium, and in-house PCR to two commercial molecular methods (the GenoType® MRSA test and the EVIGENETM MRSA Detection test) with different bacterial species in addition to S. aureus were compared. The cefoxitin disk diffusion method was superior to that of oxacillin disk diffusion and to the MIC tests in predicting mecA-mediated resistance in S. aureus when incubating at +35°C with or without the addition of NaCl to the test medium. Both the Geno Type® MRSA and EVIGENETM MRSA Detection tests are usable, accurate, cost-effective, and sufficiently fast methods for rapid MRSA confirmation from a pure culture.

Mechanisms of azole resistance and carcinogenic acetaldehyde production in chronic mucosal candidosis

Autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED, APS1) is an autoimmune disease caused by a loss-of function mutation in the autoregulator gene (AIRE). Patients with APECED suffer from chronic mucocutaneous candidosis (CMC) of the oral cavity and oesophagus often since early childhood. The patients are mainly colonized with Candida albicans and decades of exposure to antifungal agents have lead to the development of clinical and microbiological resistance in the treatment of CMC in the APECED patient population in Finland. A high incidence of oral squamous cell carcinoma is associated with oral CMC lesions in the APECED patients over the age of 25. The overall aim of this study was firstly, to investigate the effect of long-term azole exposure on the metabolism of oral C. albicans isolates from APECED patients with CMC and secondly, to analyse the specific molecular mechanisms that are responsible for these changes. The aim of the first study was to examine C. albicans strains from APECED patients and the level of cross-resistance to miconazole, the recommended topical compound for the treatment of oral candidosis. A total of 16% of the strains had decreased susceptibility to miconazole and all of these isolates had decreased susceptibility to fluconazole. Miconazole MICs also correlated with MICs to voriconazole and posaconazole. A significant positive correlation between the years of miconazole exposure and the MICs to azole antifungal agents was also found. These included azoles the patients had not been exposed to. The aim of our second study was to determine if the APECED patients are continuously colonized with the same C. albicans strains despite extensive antifungal treatment and to gain a deeper insight into the genetic changes leading to azole resistance. The strains were typed using MLST and our results confirmed that all patients were persistently colonized with the same or a genetically related strain despite antifungal treatment between isolations. No epidemic strains were found. mRNA expression was analysed by Northern blotting, protein level by western blotting, and TAC1 and ERG11 genes were sequenced. The main molecular mechanisms resulting in azole resistance were gain-of-function mutations in TAC1 leading to over expression of CDR1 and CDR2, genes linked to azole resistance. Several strains had also developed point mutations in ERG11, another gene linked to azole resistance. In the third study we used gas chromatography to test whether the level of carcinogenic acetaldehyde produced by C. albicans strains isolated from APECED patients were different from the levels produced by strains isolated from healthy controls and oral carcinoma patients. Acetaldehyde is a carcinogenic product of alcohol fermentation and metabolism in microbes associated with cancers of the upper digestive tract. In yeast, acetaldehyde is a by-product of the pyruvate bypass that converts pyruvate into acetyl-CoA during fermentation. Our results showed that strains isolated from APECED patients produced mutagenic levels of acetaldehyde in the presence of glucose (100mM, 18g/l) and the levels produced were significantly higher than those from strains isolated from controls and oral carcinoma patients. All strains in the study, however, were found to produce mutagenic levels of acetaldehyde in the presence of ethanol (11mM). The glucose and ethanol levels used in this study are equivalent to those found in food and beverages and our results highlight the role of dietary sugars and ethanol on carcinogenesis. The aims of our fourth study were to research the effect of growth conditions in the levels of acetaldehyde produced by C. albicans and to gain deeper insight into the role of different genes in the pyruvate-bypass in the production of high acetaldehyde levels. Acetaldehyde production in the presence of glucose increased by 17-fold under moderately hypoxic conditions compared to the levels produced under normoxic conditions. Under moderately hypoxic conditions acetaldehyde levels did not correlate with the expression of ADH1 and ADH2, genes catalyzing the oxidation of ethanol to acetaldehyde, or PDC11, the gene catalyzing the oxidation of pyruvate to acetaldehyde but correlated with the expression of down-stream genes ALD6 and ACS1. Our results highlight a problem where indiscriminate use of azoles may influence azole susceptibility and lead to the development of cross-resistance. Despite clinically successful treatment leading to relief of symptoms, colonization by C. albicans strains is persistent within APECED patients. Microevolution and point mutations that occur in strains may lead to the development of azole-resistant isolates and metabolic changes leading to increased production of carcinogenic acetaldehyde.

Molecular methods for the epidemiological analysis of methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae

Methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae are major health problems worldwide, both found in symptomless carriage but also causing even life-threatening infections. The aim of this thesis was to characterise MRSA and S. pneumoniae in detail by using several molecular typing methods for various epidemiological purposes: clonality analysis, epidemiological surveillance, outbreak investigation, and virulence factor analysis. The characteristics of MRSA isolates from the strain collection of the Finnish National Infectious Disease Register (NIDR) and pneumococcal isolates collected from military recruits and children with acute otitis media (AOM) were analysed using various typing techniques. Antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), spa typing, staphylococcal cassette chromosome mec (SCCmec) typing, and the detection of Panton-Valentine leukocidin (PVL) genes were performed for MRSA isolates. Pneumococcal isolates were analysed using antimicrobial susceptibility testing, serotyping, MLST, and by detecting pilus islet 1 (PI-1) and 2 (PI-2) genes. Several international community- and hospital-associated MRSA clones were recognised in Finland. The genetic diversity among MRSA FIN-4 isolates and among FIN-16 isolates was low. Overall, MRSA blood isolates from 1997 to 2006 were genetically diverse. spa typing was found to be a highly discriminatory, rapid and accurate typing method and it also qualifies as the primary typing method in countries with a long history of PFGE-based MRSA strain nomenclature. However, additional typing by another method, e.g. PFGE, is needed in certain situations to be able to provide adequate discrimination for epidemiological surveillance and outbreak investigation. An outbreak of pneumonia was associated with one pneumococcal strain among military recruits, previously healthy young men living in a crowded setting. The pneumococcal carriage rate after the outbreak was found to be exceptionally high. PI-1 genes were detected at a rather low prevalence among pneumococcal isolates from children with AOM. However, the study demonstrated that PI-1 has existed among pneumococcal isolates prior to pneumococcal conjugate vaccine and the increased antimicrobial resistance era. Moreover, PI-1 was found to associate with the serotype rather than the genotype. This study adds to our understanding of the molecular epidemiology of MRSA strains in Finland and the importance of an appropriate genotyping method to be able to perform high-level laboratory-based surveillance of MRSA. Epidemiological and molecular analyses of S. pneumoniae add to our knowledge of the characteristics of pneumococcal strains in Finland.